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ABSTRACT 

Self-Interacting Dark Matter is an attractive alternative to the Cold Dark Matter paradigm only 
if it is able to substantially reduce the central densities of dwarf-size haloes while keeping the 
densities and shapes of cluster-size haloes within current constraints. Given the seemingly 
stringent nature of the latter, it was thought for nearly a decade that Self-Interacting Dark 
Matter would be viable only if the cross section for self-scattering was strongly velocity- 
dependent. However, it has recently been suggested that a constant cross section per unit mass 
of cry/m ~ 0. 1 cm 2 g~ 1 is sufficient to accomplish the desired effect. We explicitly investigate 
this claim using high resolution cosmological simulations of a Milky-Way size halo and find 
that, similarly to the Cold Dark Matter case, such cross section produces a population of 
massive subhaloes that is inconsistent with the kinematics of the classical dwarf spheroidals, 
in particular with the inferred slopes of the mass profiles of Fornax and Sculptor. This problem 
is resolved if ar/m^ 1 cm 2 g _1 at the dwarf spheroidal scales. Since this value is likely 
inconsistent with the halo shapes of several clusters, our results leave only a small window 
open for a velocity-independent Self-Interacting Dark Matter model to work as a distinct 
alternative to Cold Dark Matter. 
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1 INTRODUCTION 

It is now clear that observations of dark matter dominated sys- 
tems such as low-mass (dwarf) and low surface brightness (LSB) 
galaxies favour the presence of dark matter cores of 0(1 kpc) (e.g. 
|Moore|1994||Kuzio de Naray et al.|2008||de Blok|2010||Walker&l 
|Penarrubia|201 l||Amorisco* & Evans 2 012) . These observations are 
a challenge for the Cold Dark Matter (CDM) paradigm where dark 
matter haloes are predicted to have density cusps, an imprint of the 
collisionless nature of CDM (the core-cusp problem). In a possibly 
related issue, it has been pointed out recently that the dark satellites 
of Milky-Way (MW) size halo simulations are too dense to be con- 
sistent with the kinematics of the MW dwarf spheroidals (dSphs) 
(the too big to fail problem; Boylan-Ko lchin et al.|201 1[|2012) . 

Although these are significant challenges to the CDM model, 
their solution could naturally lie in our incomplete understanding of 
the complex process of galaxy formation. Even though the internal 
dynamics of dwarfs is dominated by dark matter today, it is con- 
ceivable that earlier episodes of star formation and subsequent gas 
removal by supernova feedback might have been violent enough to 
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modify the initial cuspy dark matter distribution into a cored one 
(e.g. Navarro et al. 1996, Pontzen & Govemato 2012). Current hy- 
drodynamical simulations have shown that such a mechanism is 
able to create large cores in intermediate-mass galaxies ( |Governato| 
|et al.|2 010 2012), and, once tidal snipping is taken into account, 
it might solve the too big to fail problem as well (Brooks & Zolo- 
|tov||2012] >. It is not clear, however, if such episodes of large gas 
blow-outs are consistent with the star formation histories and stel- 
lar properties of LSBs and dwarfs (e.g. Kuzio de Naray & Spekkens 



[20fr]|Boylan-Kolchin et al.|2012||Penarrubia et al.|2012) 

Given the large uncertainties regarding whether baryon 
physics can reconcile the CDM model with observations of dwarf 
galaxies, it is prudent to consider the alternative, which is to ques- 
tion the fundamental CDM hypotheses, namely, the collisionless 
and cold nature of CDM particles. This alternative is additionally 
encouraged by the null detection of several experiments that are 
pursuing the discovery of the favoured CDM particles, and whose 
sensitivity is reaching the natural values for the interaction cross 
sections of the particle physics models that predict them (e.g. Su- 
persymmetry, Abazajian et al. 2012, A prile et al.|2012) . 

An exciting possibility is that of self-interacting dark mat- 
ter (SIDM) originally introduced over a decade ago by Sp ergel &| 
Steinhardt (2000). Self-scattering between dark matter particles is 
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a feature of present hidden-sector dark matter models that predict 
the existence of new gauge bosons. The presence of these bosons is 
invoked to enhance the annihilation and/or self-scattering of dark 
matter particles to explain a number of puzzling observations (the 
Sommerfeld enhancement, e.g. |Arkani-Hamed et al. 2009 Buck-| 
ley & Fox 2010). Collisional dark matter is constrained by the re- 
quirements from different astrophysical observations, such as the 
ellipsoidal shape of haloes, the avoidance of subhalo evaporation 
in galaxy clusters, and the avoidance of the gravothermal catastro- 
phe (e.g.|Miralda-Escu de 2002, Gnedin & Ostri ker|200 1 1 |Firmani| 
|et al.|200l| . The original excitement caused by SIDM died off by 
the apparently strong constraints on the scattering cross section set 
particularly by X-ray and lensing observations of clusters in the 



analysis by Miralda-Escude (2002): or/msC 0.02 cm g ; such 



a low cross section would have no relevant impact for the dynam- 
ics of galaxies at the 0(1 kpc) scale. [Peter et al.| ( |2012| > have re- 
vised this constraint and found that it was overestimated by over 
an order of magnitude suggesting that a current constraint is of 
O(0.1 cm 2 g- 1 ). 

SIDM is clearly a viable model if the cross section depends 
on the relative velocity in such a way that dark matter behaves as 
a collisional fluid in dwarfs, and is essentially collisionless at the 
scale of clusters. Although this idea was phenomenologically pro- 
posed ( Yoshida et al. 2000 1 and explored with cosmological simu- 
lations l Colin et al.|2002) a decade ago, its theoretical support has 



come up only recently (e.g. 


Ackermanetal.|2009, 


Feng et al. 2009 


2010||Buckley & Fox|2010 


|Loeb&Weiner|2011| 


van den Aarssen 


et al. 2012, Tulin et al. 2012 ). Moreover, earlier simulations lacked 



the resolution needed to reliably explore the sub-kpc region of 
dwarf-size haloes. It was only until recently that it was explicitly 
shown that theoretically motivated velocity-dependent SIDM (vd- 
SIDM) models produce core sizes consistent with those found in 
MW dSphs, and also solve the emergent too big to fail problem 



Vogelsberger et al.|2012 


hereafter VZL). 


On the other hand, 


Rocha et al.l (12012b have suggested that 



a velocity-dependent cross section is not essential since a SIDM 
model with a constant ut /m = 0.1 cm 2 g _1 (allowed by cluster 
constraints) is also consistent with the inner structure of dSphs. 
This seems to contradict earlier estimates made by Yoshi da et al.| 
(2000) who suggested that the average number of collisions per 
particle in the central core, n corc , scales as the cube root of the halo 
mass. Since for this cross section, n corc ~2 for a cluster-size halo 
(within ~ 10% of the virial radius, see their Fig. 2), n corc would 
thus be suppressed by a factor of 0(100) in dwarf-size haloes, re- 
sulting in cores that are too small. This scaling is however imprecise 
since for a halo of virial mass M: 



'(pcorc) [ — ) (o"vcl)ta 

\ m I 



V m I 



/3, 
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where (p C ore) and (cr ve i) are the average density and velocity dis- 
persion of dark matter particles, t agc is the formation time, and 
A coro is proportional to the density contrast relative to the back- 
ground density; all of these quantities are defined within the core of 
the halo. Thus, the higher concentrations and larger formation times 
of dwarf-size haloes reduce the M 1 ^ 3 dependence. It is important 
to remark that both, |Yoshida et al.| ( |2000l and |Rocha et al.| ( |2012) , 
extrapolated the regimes they could directly simulate to the regime 
of dwarfs. In this Letter, we resolve this issue by proving explic- 
itly that a constant scattering cross-section of ot /m = Q.\ cm 2 g _1 
is not able to create 0(1 kpc) cores in the dark subhaloes where 
the MW dSphs are expected to live; it deviates only slightly from 
the CDM predictions. Unless baryonic processes are invoked, the 
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Figure 1. Dependence of the momentum-transfer weighted cross-section 
per unit mass on the relative velocity for the different SIDM models consid- 
ered here. The constant cross section cases with o^/m^ 1 cm 2 g _1 are 
likely ruled out by halo shapes based on X-ray and lensing observations of 
clusters (Peter et al. 20121. The models with a velocity-dependent cross- 
section are tuned to satisfy all current astrophysical constraints and have 
been shown to be consistent with the kinematics of MW dSphs (VZL). 



range of interesting constant a^/m values is thus very narrow: 
0.1 cm 2 g -1 <<tt /m< 1 cm 2 g _1 . 



2 SIMULATIONS AND RESULTS 

Our analysis is based on re-simulations of the Aq-A halo (level 3 
resolution) of the Aquarius project ( Springel et al. 2008 1, which is a 
set of representative MW-like haloes within the CDM WMAP-lyr 
cosmology. This halo has a virial mass of Af 2a o~1.8 x 1O 12 M 
within a radius of 7-200 ~ 246 kpc (enclosing an average density of 
200 times the critical density). The particle mass in the simulations 
is m p ~4.9 x 10 4 M and the Plummer equivalent gravitational 
softening length is e~ 120 pc. We use an algorithm that adds dark 
matter self-scattering to the A^-body code GADGET-3 for gravita- 
tional interactions (last described in Springel 2005). The algorithm 
uses a N-body/Monte Carlo approach to represent the microphysi- 
cal scattering process in the macroscopic context of the simulation. 
The connection between this type of approach and the Boltzmann 
equation is nicely described in Appendix A of R ocha et al"| ( |2012] >. 
The details of this algorithm can be found in VZL, as well as sim- 
ple controlled tests that show the agreement between the outcome 
of the code and analytical expectations. All dark matter models 
were simulated starting with the same initial conditions and their 
present-day self-bound subhalo population was identified using the 
SUB FIND algorithm JSpringel et al.|2001) . 

In addition to CDM, we consider five SIDM cases (recently 
presented in Vogelsberger & Zavala 2012 to analyse the impact of 
self-scattering in direct detection experiments): three with a con- 
stant cross section and two with a velocity-dependent one given by 
a Yukawa-like interaction (e.g. Loeb & Weiner 201 1 1. The transfer 
cross section scaling with the relative velocity can be seen in Fig. [I] 
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Figure 2. The circular velocity profiles at z = encompassing the 1 st and 
3 rd quartiles of the distribution of the 15 subhaloes with the largest values 
of V ma x(z = 0). The symbols with error bars are estimates of the circu- 
lar velocity within the half-light radii for 9 MW dSphs (Wal ker et al.|2009[ 
|Wolf et al.|2010) . Clearly, the most massive CDM subhaloes are inconsis- 
tent with the kinematics of the MW dSphs. SIDM can alleviate this problem 
only for a constant scattering cross-section /m > 1 cm 2 g _1 (SIDM 10 
and SIDM1) or if it has a velocity dependence (vdSIDMa and vdSIDMb). 
Current constraints from clusters put an upper limit to the constant cross 
section case close to ar/m~0.1cm 2 g — 1 (SIDM0.1). This value is too 
low to solve the too big to fail problem. The observational data in the bot- 
tom right can be fitted by lower mass subhaloes, not shown here since they 
are affected by the limited resolution of our simulations. 



We note that the formula for ar/m for the velocity-dependent 
cases is only valid in the classical regime, once quantum effects 
are important, the finite interaction length of the Yukawa poten- 
tial cuts off the zero-velocity divergence of the cross section (see 
e.g. |Feng et al.|20T0| >. For our purposes, the quantity of relevance 
is (ctt/w.) v which goes to zero at zero velocity. It is clear that 
for the vdSIDM models, UT/m> 0.1 cm 2 g _1 at the character- 
istic velocities in MW dSphs (the observed v elocity dispersion o f 
stars along the line of sight is ~ 10 km s _1 , e.g. |Walker et al.|2009| . 
This fact alone already casts a doubt on the possibility of SIDM0. 1 
(or/m = 0.1 cm 2 g _1 ) producing similar results as the vdSIDM 
cases that were shown to be consistent with the kinematics of the 
MW dSphs in VZL. We note that there is a change in nomen- 
clature relative to VZL: RefP0=CDM, RefPl=SIDM10, RefP2- 
3=vdSIDMa-b. 

Fig|2] shows the inter-quartile range (i.e., 25-75%) of the dis- 
tribution of the present-day circular velocity profiles of the 15 sub- 
haloes with the largest values of Vmax(z = 0) (the maximum 
of the circular velocity) within 300 kpc halocentric distance. The 
symbols with error bars correspond to estimates of the circular 
velocity within the half-light radii of the sample of 9 MW dSphs 
used by Boylan-Kolchin et al. (2011, 2012). Since current data for 
the stars in the dSphs provide an incomplete description of the 6- 
dimensional phase-space distribution, the derived mass profiles are 
typically degenerate with the velocity anisotropy profile. However, 



the uncertainty in mass that is due to this degeneracy is minimised 
near the half-light radius, where Jeans models tend to give the same 
value of enclosed mass regardless of anisotropy (e.g. |Strigari et al.| 
[20071 |Walker et al.|2009||Wolf et al.|20T0> . Observations can then 
be used to constrain the maximum dark matter density within this 
radius. CDM clearly predicts a population of massive subhaloes 
that is inconsistent with all the 9 dSphs, whereas for SIDM this 
problem disappears as long as (JT/m> 1 cm 2 g _1 on dSph scales. 
The currently allowed case with <tt /m = 0.1 cm 2 g _1 is very close 
to CDM, only reducing slightly the inner part of the subhalo veloc- 
ity profiles. On the contrary, the vdSIDM models clearly solve the 
too big to fail problem. We note that the extent of the too big to fail 
problem in CDM depends on the mass of the MW halo, if it is in 
the low end of current estimates, < 10 12 Mq, the problem may be 
resolved (e.g. Wang et al. 2012 ), although a low halo mass may gen- 
erate other difficulties such as explaining the presence of the Mag- 
ellanic Clouds. In the context of SIDM, the lower the mass of the 
MW halo, the weaker the argument against ar /m = 0.1 cm 2 g _1 . 

A simple statistical test of the agreement between the subhalo 
distributions of two models and the 9 dSphs is to compute the chi- 
square difference associated to the likelihood of having n + (n~) 
data points above (below) the median of the distribution of each 
model. Assuming that the probability distribution of finding n 
data points is Poissonian: 

A X 2 = 2(ln(n+! n^l) - ln(n+! nj!)). (2) 

Comparing SIDM1 and the vdSIDM models with SIDM0.1, 
the difference is driven solely by Draco with the former pre- 
ferred over the latter with A% 2 ~4.4 (2.1<r). Using an interpo- 
lation of our three constant cross section cases, we estimate that 
ar/m^ 0.6 cm 2 g -1 is the minimum value for which A\ 2 = 
relative to SIDM1. 

To show the typical core size and central densities that are pre- 
dicted by the different SIDM models, we plot in Fig.[3]fhe density 
profile of the 15 subhaloes with the largest Vm ax (z = 0) values. 
A value of ar/m^l cm 2 g _1 is needed for a constant cross sec- 
tion SIDM model to mimic the effect of the vdSIDM models and 
produce ~ 1 kpc cores with central densities of O(0.1 Mq pc~ 3 ). 
If the transfer cross section is reduced to 0.1cm 2 g _1 , then the 
subhaloes are only slightly less dense than in CDM, having cores 
(central densities) that are at least twice smaller (higher) than those 
in the other SIDM cases. 

VZL showed that the SIDM 10 and vdSIDM models have con- 
vergent density and circular velocity profiles within the central den- 
sity core; we have found the same for SIDM1 and to lesser extent 
for SIDM0.1. Convergence is harder to achieve for CDM since, at a 
fixed radius, the two-body relaxation time is shorter than for SIDM 
(due to the reduced densities in the latter case). [Power et al.|(2003| > 
showed that the density profile converges at a given radius when 
the two-body relaxation time is larger than the Hubble time at this 
radius. At the resolution level of our simulations, the convergence 
radius for CDM is ~ 600 pc, which implies that the CDM circular 
velocity and density profiles shown in Figs.[2]and[3]underestimate 
the true dark matter content within ~ 600 pc ( Springel et al. 2008 1, 
whereas for SIDM is at least half of this value. In any case, the ex- 
pectation is that if the density profile of SIDM0. 1 has not converged 
yet, higher resolution would drive it towards higher densities, not 
lower, bringing it even closer to CDM (this is a trend confirmed for 
the cases analysed in VZL, see their Fig. 9). 

By using the fact that some MW dSphs have chemo- 
dynamically distinct stellar subcomponents that independently 
trace the same gravitational potential, Walker & Penarrubia {201 1| | 
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Figure 3. Density profile of the 15 subhaloes with the largest V m3x {z = 
0) values within CDM and different SIDM models (see Fig. [TJ. We 
show the median and 1 st and 3 rd quartiles of the subhalo distribution 
for each case. The velocity-dependent SIDM cases produce cores of ap- 
proximately 600 pc. Of the constant cross section SIDM models we ex- 
plored, the one that is currently allowed by cluster constraints, SIDM0.1 
(<tt /m = 0.1 cm 2 g _1 ), only deviates slightly from CDM; the associated 
core sizes are less than 300 pc. 

showed that it is possible to constrain the slopes of their inner mass 
profiles. They found that Fornax and Sculptor are consistent with 
cored density profiles while cuspy profiles with p oc r~ are ruled 
out with a significance > 96% and > 99%, respectively. We use 
this method to test the consistency of the SIDM models explored 
here. We found that all SIDM models, except for SIDM0.1, are 
well fitted by the following three-parameter formula: 

(r + r c ) (r 2 + rj) 

which is similar to the Burkert profile (Burkert 1995 1 but with two 
scale radii r s and r c . The remaining case, SIDM0.1, is better fitted 
by: 

I \ = Po r 3 s 

P(T) ~ (r + r c )(r + r s ) 2 - (> 

Using these formulae, we found the best fit parameters for the mas- 
sive subhaloes in each of the SIDM models. Such fits are restricted 
to a radial range between the softening length of our simulations 
~ 120 pc and the radius where tidal stripping has made the outer 
logarithmic slope of the density profile steeper than —3. The latter 
restriction is of relevance only for four subhaloes that are affected 
significantly by tidal stripping within ~ 5 kpc. Two of these are 
clearly affected within ~ 1 kpc and should likely be removed in a 
more detailed analysis; they are the least consistent with the data. 
To find the best fit parameters we minimise: 

Q 2 = ^— E Q*Pi(ri) - lnpflUn)) 2 , (5) 

% 

where the sum goes over all radial bins.Thus defined, Q gives an 
estimate of the goodness of the fit. In Table [T] we give the best 



Name p(r = 200 pc) [Mq kpc 3 ] r s [kpc] r c [kpc] 



vdSIDMa 


1.37 x 10 s 


0.94 


0.75 


vdSIDMb 


1.37 x 10 8 


0.94 


0.73 


SIDM1 


1.16 x 10 8 


0.96 


1.33 


SIDM0.1 


2.31 x 10 8 


0.97 


0.41 



Table 1. Best fit parameters for the median of the SIDM density profiles of 
the 15 subhaloes with the largest V ma x(z = 0) values. The last two have 
a constant cross section while the others have a velocity-dependent cross 
section (see Fig.fTJ. SIDM1.0 is likely ruled out by cluster observations 
(see |Rocha et al.|2012) . The density profile used for the fits is given by 
Eq. {3} for SIDM1 and the vdSIDM models, and by Eq. (3) for SIDM0.1. 

fit parameters for the median of the subhalo population for each 
SIDM model (except for SIDM10 which has been ruled out). We 
note that Penarr ubia et al.| ([2012l already used Eq. {4} to estimate 
r c > 1 kpc for Fornax and Sculptor. Cores of this size are too large 
to be consistent with most of the subhaloes in SIDM0. 1 . 

To test the consistency of the different SIDM models, we 
use the parameters of the fits to compute the slope of the inner 
mass profile between the pair of half light radii (the median like- 
lihood values) of the two distinct stellar subcomponents in For- 
nax and Sculptor. We then test whether this slope is as steep as 
the lower limit set by the data. The confidence level at which 
a given slope is said to be excluded is determined by the frac- 
tion of the posterior distribution f p of allowed slopes that are 
larger. For a? /m = 0.1 cm 2 g _1 , all but 2 subhaloes are excluded 
at > 95(90)% confidence for Fornax (Sculptor); the remaining two 
subhaloes have values of f p > 0.86(0.81) for Fornax (Sculptor). 
On the contrary, the other SIDM models (except for SIDM 10 that 
was not analysed) are clearly more consistent with the data with 
only four subhaloes excluded at 90% confidence for Fomax (five of 
the subhaloes actually have f p < 0.8), while only three subhaloes 
are excluded at > 80% confidence for Sculptor. We found no clear 
preference between the vdSIDM models and the case with constant 
Ot /m — 1 cm 2 g _1 . To consider the impact of the non-spherical 
morphologies of Fornax and Sculptor, we repeated the analysis for 
elliptical rather than circular radii for the stars used to estimate the 
slope of the mass profiles (see sect. 6.1 of Walker & Penarrubia| 
|201 1[ >. We find that for all models Fornax becomes slightly more 
exclusive while Sculptor is considerably less exclusive. 



3 DISCUSSION AND CONCLUSIONS 

Self-Interacting Dark Matter (SIDM) offers a promising solution 
to the dwarf-scale challenges faced by the otherwise-remarkably 
successful Cold Dark Matter (CDM) model. The original idea of a 
velocity-independent, elastically scattering cross section died off 
quickly, mostly due to the apparently stringent constraint found 
by pira lda-Escude"| {2002} requiring that the cross-section per unit 
mass was ot jm ^ 0.02 cm 2 g~ 1 . This value is uninteresting, with 
earlier estimates requiring or/m to be at least of 0(1 cm 2 g _1 ) 
to create ~ 1 kpc cores in dwarf-size haloes ( |Yoshida et al.|2000"| 
|Dave et al.|2001[ >. |Peter et"a l. (2012) have recently revised earlier 
constraints on collisional dark matter and found them to be over- 
estimated by over an order of magnitude; the current constraint 
is <jt /m<0.1 cm 2 g _1 . Moreover, these authors have revived, in 
a companion paper ( |Rocha e t al. 2012}, the velocity-independent 
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Figure 4. Subhalo mass function for a MW-size halo within CDM and dif- 
ferent elastic SIDM models. The only model that leads to a difference rel- 
ative to CDM has a constant cross section of oj-/m = 10 cm 2 g _1 , which 
is clearly ruled out by cluster observations. 

SIDM model by suggesting that a value of ot /m = 0.1 cm 2 g _1 is 
seemingly consistent with the inner structure of the MW dSphs. 

Motivated by the prospect of a viable constant cross section 
SIDM model, we investigate the claims from |Rocha et al.|j2012) 
using high resolution cosmological SIDM simulations of a MW- 
size halo. Contrary to Rocha et al. (2012), we are able to resolve 
the sub-kpc structure of the massive subhalo population to suffi- 
ciently small radii for comparison with the MW dSphs. We find that 
a velocity-independent SIDM model is consistent with the kinemat- 
ics of dSphs only if ot /m~ 1 cm 2 g _1 (see Fig.[2]l, i.e., a value of 
this order is required to solve the too big to fail problem (Boylan- 
|Kolchin et al. 201 1 , 2012). If the cross section is lower by an order 
of magnitude, the subhalo population is still too dense to be con- 
sistent with the MW dSphs. On the other hand, as shown already 
in VZL, velocity-dependent SIDM models with a Yukawa-like in- 
teraction (as proposed in |Loeb &" Weiner 201 1 see Fig.^l success- 
fully solve the too big to fail problem. 

We also use the inner slopes of the mass profiles of For- 
nax and Sculptor, from the analysis of Walker & Pen arrub?a| 
< |20 1 1 1 >, as examples to test the consistency of the different mod- 
els we simulate here. For a velocity-independent SIDM model with 
ot /m~0.1 cm 2 g" 1 , we find that 13 of of the 15 subhaloes with 
the largest V ma:K (z — 0) values are inconsistent with the data from 
Fornax (Sculptor) at > 95(90)% confidence (the other two are in- 
consistent at > 81% confidence). A constant cross section ten times 
larger is as consistent as the velocity-dependent SIDM models ex- 
plored here with only four (three) of the top 15 subhaloes excluded 
at > 90(80)% confidence in the case of Fornax (Sculptor); for 
all these cases, there are several subhaloes that are unambiguously 
consistent with the data. 

According to the analysis of |Peter et al.| (2012 >, a constant 
cross section of ot /m — 1 cm 2 g~ is likely inconsistent with the 
observed halo shapes of several clusters. We have now shown that 
o~t /m = 0.1 cm 2 g^ 1 is too close to CDM to represent a distinct 
alternative. An interpolation of our simulations suggests that the 
central densities of the massive subhaloes would be consistent with 
the MW dSphs if <7T/m~0.6cm 2 g _1 . We conclude that the hy- 
pothesis of a constant scattering cross section as solution to the 
core-cusp problem remains viable but within a very narrow range 



of ot /m values. The challenges to make a definitive test of this hy- 
pothesis are twofold: the cluster-constraints need be refined, and the 
impact of conservative baryonic processes needs to be estimated. 
Although adding gas physics is the next step of SIDM simulations, 
a challenge to make SIDM an even more attractive alternative to 
CDM is the prospect of explaining the observed scarcity of MW 
satellites and field dwarfs (e.g. Klypin et al. 1999; Zaval a~et al.| 
2009 ) without invoking extreme baryonic processes. As we show 
in Fig. [4] all allowed elastic SIDM models essentially produce the 
same abundance of dwarf-size haloes as in CDM. A promising pos- 
sibility is that of exothermic interactions between excited and non- 
excited states of dark matter (e.g. |Loeb & W einer 2011). The ve- 
locity kick imparted during the collision might be large enough to 
cause the evaporation of low-mass haloes. 
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